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HIGHLIGHTS 


•  A  transient  TEG  model  was  proposed  with  consideration  of  heat-electricity  coupling. 

•  The  dynamic  characteristics  of  the  TEG  was  studied  numerically. 

•  Four  kinds  of  variations  in  application  environments  were  designed. 

•  Response  hysteresis  of  power  and  overshoot/undershoot  of  efficiency  were  found. 
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The  practical  application  environments  of  thermoelectric  generators  (TEGs)  always  change,  which 
make  a  requirement  for  studying  the  dynamic  response  characteristics  of  TEGs.  This  work  develops  a 
complete,  three-dimensional  and  transient  model  to  investigate  this  issue.  The  model  couples  the  energy 
and  electric  potential  equations.  Seebeck  effect,  Peltier  effect,  Thomson  effect,  Joule  heating  and  Fourier 
heat  conduction  are  taken  into  account  in  this  model.  Dynamic  output  power  and  conversion  efficiency 
of  the  TEG,  which  are  caused  by  variations  of  the  hot  end  temperature,  cold  end  temperature  and  load 
current,  are  studied.  The  response  hysteresis  of  the  output  power  to  the  hot  end  and  cold  end  temper¬ 
atures,  the  overshoot  or  undershoot  of  the  conversion  efficiency  are  found  and  attributed  to  the  delay  of 
thermal  diffusion.  However,  the  output  power  is  synchronous  with  the  load  current  due  to  much  faster 
electric  response  than  thermal  response. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Alternative  green  energy  technology  which  aims  to  reduce 
greenhouse  gas  emissions  and  fossil  fuel  usage  has  become  one  of 
the  most  popular  researches  in  the  background  of  the  current  en¬ 
ergy  crisis.  As  one  of  the  alternatives,  thermoelectric  technology 
has  attracted  great  attention  in  the  last  decade  due  to  its  excellent 
features,  such  as  no  mechanical  moving  parts,  no  refrigerants, 
layout  flexibility,  low  maintenance  costs,  long  run-time  and 
other  advantages  [1—3].  Thermoelectric  devices,  which  consist  of 
tens  to  hundreds  of  pairs  of  p-type  and  n-type  semiconductors,  can 
be  divided  into  two  different  groups:  thermoelectric  cooler  (TEC) 
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[4-6],  which  directly  convert  electricity  to  thermal  energy  for 
cooling  purpose  based  on  the  Peltier  effect,  and  thermoelectric 
generator  (TEG)  [7-9],  driven  by  large  temperature  differences  to 
convert  the  heat  into  electrical  power  based  on  the  Seebeck  effect. 

The  TEG  performance,  including  the  output  power  and  conver¬ 
sion  efficiency,  is  not  only  related  to  the  properties  of  thermo¬ 
electric  materials,  but  also  closely  related  to  the  operating 
conditions.  For  a  long  time,  the  TEG  conversion  efficiency  had  been 
very  low  due  to  limitation  of  thermoelectric  materials,  hence,  the 
TEG  technology  was  only  used  in  aerospace  and  other  cutting-edge 
areas  [10].  However,  with  the  successful  development  of  new 
thermoelectric  materials  and  the  significant  improvement  of  the 
TEG  conversion  efficiency,  the  TEG  technology  has  been  regarded  as 
a  promising  power  generator  in  civilian  applications,  thus,  a 
number  of  works  have  been  carried  out  to  investigate  the  charac¬ 
teristics  of  the  TEG  [11-20]. 

Almost  all  of  the  existing  literature  of  TEGs  carried  out  the 
performance  evaluation  in  steady-state  conditions  [7-9,11-20], 
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the  transient  performance  was  concerned  only  in  a  few  litera¬ 
ture  [21-23].  In  fact,  the  practical  application  environments  of 
TEGs  always  change.  First  of  all,  the  heat  source  temperature  will 
change  with  varied  supplies,  for  example,  the  coal  or  fuel  supply 
rate  can  not  be  static  in  iron  and  steel  works,  power  generating  and 
chemical  industries,  leading  to  time-dependent  industrial  waste 
heat;  for  the  collection  of  automobile  exhaust  heat,  the  amount  of 
heat  will  increase  with  the  acceleration  and  reduce  with  the 
deceleration  of  the  automobile.  Secondly,  the  consumption  de¬ 
mand  for  electrical  energy  inevitably  varies  in  different  time  scale, 
which  makes  a  requirement  for  the  change  of  load  current  in  the 
TEG  system.  Thirdly,  the  temperature  of  the  cold  end  of  the  TEG  will 
be  changed  with  the  flow  rate  of  cooling  medium  in  heat  sinks. 
Furthermore,  the  thermoelectric  transient  behavior  is  particularly 
apparent  in  the  start-up  and  shut-down  stages  of  the  devices.  Once 
the  above-mentioned  changes  occur,  the  output  power  and  con¬ 
version  efficiency  of  the  TEG  will  change,  showing  the  dynamic 
behavior,  the  original  steady-state  equilibrium  is  destroyed,  while 
the  formation  of  a  new  steady-state  will  inevitably  consume  a 
certain  time. 

Chen  et  al.  [21  ]  studied  the  transient  behavior  of  the  TEG  when 
the  user  load  was  changed  based  on  the  energy  balance  equations 
only  at  the  hot  and  cold  ends  (zero-dimensional  model),  the  model 
divided  both  the  Joule  heat  and  the  Thomson  heat  into  two  equal 
parts,  one  part  flows  to  the  hot  end,  and  the  other  to  the  cold  end  of 
the  TEG,  however,  it  is  not  strictly  accurate,  because  the  Joule  heat 
and  Thomson  heat  are  the  inner  heat  sources.  Moreover,  this 
method  does  not  solve  the  internal  temperature  distribution  of  the 
TEG,  so  that  dependence  of  TEG  performance  on  the  internal  heat 
transfer  and  carrier  transport  can  not  be  understood  well.  Mon- 
tecucco  et  al.  [22]  used  the  method  of  solving  one-dimensional  heat 
conduction  equation  to  study  the  transient  characteristics  of  the 
TEG  when  the  load  was  changed,  but  ignored  the  three- 
dimensional  effect  of  the  TEG.  Alata  et  al.  [23]  also  developed  a 
one-dimensional  transient  TEG  model  to  investigate  the  dynamic 
variation  of  conversion  efficiency  with  increased  load  current. 
Compared  with  Montecucco  et  al.’s  model,  the  Thomson  effect  was 
taken  into  account  in  Ref.  [23].  In  these  models  [21-23],  the  current 
density,  J,  was  assumed  to  be  distributed  uniformly  within  the 
whole  TEG,  so  that  J  =  //A,  where  I  is  the  load  current  and  A  is  the 
cross-sectional  area  of  semiconductors.  However,  our  previous 
[24,25]  and  Cheng’s  [26]  works  have  shown  that  there  exists 
coupling  of  temperature  field  and  electric  potential  field  within 
semiconductors  for  the  TEC,  which  causes  three-dimensional 
temperature  and  current  distributions,  especially  for  high  tem¬ 
perature  differences  and  large  applied  currents.  Thus,  it  is  expected 
that  this  coupled  effect  may  also  be  important  for  the  TEG,  because 
the  TEG  generally  operates  in  larger  temperature  differences  than 
the  TEC  in  order  to  produce  a  higher  output  power.  Furthermore, 
these  studies  [21-23]  focused  on  the  dynamic  characteristics  only 
when  the  load  current  was  changed,  in  practice,  transient  operating 
conditions  of  the  TEG  are  caused  by  variation  of  not  only  the  load 
current,  but  also  temperatures  of  the  hot  or  cold  ends.  Thus,  un¬ 
derstanding  the  dynamic  characteristics  of  the  TEG  induced  by 
variations  of  hot  or  cold  heat  sources  is  necessary  for  the  design  and 
operation  of  the  TEG. 

The  objective  of  this  work  is  to  develop  a  fully  three- 
dimensional  and  transient  TEG  model,  which  introduces  the 
coupling  of  heat  transfer  and  electric  conduction  within  semi¬ 
conductors.  The  model  is  adopted  to  investigate  the  transient 
characteristics  of  the  TEG.  Four  kinds  of  variation  patterns  (i.e.  step 
increase,  linear  increase,  step  decrease  and  linear  decrease)  in  the 
hot  end  temperature,  the  cold  end  temperature  and  the  load 
current  are  respectively  designed  to  reveal  the  dynamic  output 
power  and  conversion  efficiency  of  the  TEG,  which  can  provide 


some  useful  insight  into  response  mechanism  for  practical  TEG 
applications. 

2.  Model 

A  TEG  device  in  essence  is  a  thermopile  composed  of  a  number 
of  p-  and  n-type  semiconductor  pairs.  Because  the  thermal  char¬ 
acteristics  are  the  same  for  all  TEG  elements,  only  one  element  is 
considered  here.  Fig.  1  shows  the  schematic  of  the  TEG  element,  it 
includes  a  columnar  p-type  semiconductor  and  a  columnar  n-type 
semiconductor  which  are  sandwiched  between  copper  connectors, 
there  are  two  ceramic  plates  at  both  ends  of  the  connectors  for 
insulation.  The  thicknesses  of  the  ceramic  plate,  metal  connector 
and  semiconductor  are  Ho,  Hi  and  H2,  respectively,  the  cross- 
sectional  areas  of  p-  and  n-type  semiconductors  are  same  and 
equal  to  I2  x  I2  with  a  separation  distance  of  Li. 

Charge  carriers  in  p-  and  n-type  semiconductor  are  holes  and 
electrons,  respectively.  When  the  TEG  element  is  sandwiched  be¬ 
tween  the  hot  and  cold  sources,  the  hot  end  of  the  TEG  absorbs  heat 
from  the  hot  source  and  the  cold  end  releases  heat  to  the  cold 
source  due  to  carriers  transport,  which  causes  conversion  of  ther¬ 
mal  energy  into  electrical  energy  in  the  form  of  electromotive  force 
or  output  current. 


2.1.  Governing  equations 


The  three-dimensional  transient  TEG  model  established  in  this 
paper  includes  the  energy  equation  and  the  electric  potential 
equation.  The  energy  equations  of  the  metal  connector,  the  p-type 
semiconductor,  the  n-type  semiconductor,  and  the  ceramic  plate 
are  as  follow: 


V-  (A/VT) 


fti-VT 


(1) 


where,  p  is  the  density,  cp  is  the  specific  heat,  A  is  the  thermal 
conductivity,  a  is  the  electric  conductivity,  and  0  is  the  Thomson 
coefficient.  The  subscript  i  is  conn  for  connectors,  p  for  the  p-type 
semiconductor,  n  for  the  n-type  semiconductor,  cer  for  ceramic 
plates,  respectively.  The  first  term  on  the  right  side  in  Eq.  (1)  de¬ 
notes  the  Fourier  heat  conduction,  the  second  and  third  terms 
denote  the  internal  heat  sources  due  to  the  Joule  heating  and 
Thomson  effect  J(x,y,  z)  is  the  local  current  density,  and  J  =  0  for  the 
ceramic  plate.  The  Thomson  coefficient  is  proportional  to  the  first 
derivative  of  the  Seebeck  coefficient  versus  temperature,  or: 


- 

d T 


(2) 


where,  a  is  the  Seebeck  coefficient  of  semiconductors. 

The  driving  force  for  carriers  is  the  electric  potential,  which  can 
be  expressed  as  follows: 


V- 


(a(V(j>  -  aWT)  +  eV 


9 1) 


=  0 


(3) 


where,  0  is  the  electric  potential,  aiVT  is  the  Seebeck  electromotive 
force  coming  from  the  Seebeck  effect,  e  is  the  permittivity. 

The  electric  field  can  be  calculated  by  the  following  equation: 

t  =  -V0  +  aVT  (4) 

Thus,  the  current  density  in  Eq.  (1 )  can  be  calculated  as  follows: 


J  =  at 


(5) 
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Generally,  the  performance  of  the  TEG  is  evaluated  by  two  pa¬ 
rameters:  output  power,  P,  and  conversion  efficiency,  rjf  which  are 
defined  as  follows: 


P  =  IV 

(6) 

p 

(7) 

where,  /  is  the  load  current,  V  is  the  output  voltage,  Qh  is  the  heat 
supplied  to  the  hot  end  of  the  TEG. 

2.2.  Initial  and  boundary  conditions 

The  purpose  of  this  paper  is  to  study  the  dynamic  characteris¬ 
tics  of  the  TEG  when  the  application  environments  are  changed, 
including  hot  end  temperature,  cold  end  temperature  and 
load  current.  The  boundary  conditions  of  the  model  are  set  as 
follows. 


For  a  practical  operation,  the  TEG  absorbs  heat  from  a  hot  source 
such  as  industrial  flue  gas  or  automobile  exhaust  gas,  and  the  heat 
is  dissipated  by  a  heat  sink  (cold  source).  Thus,  the  third  kind 
boundary  condition  should  be  used  if  the  hot  and  cold  sources  are 
not  included  in  the  modeling.  However,  for  simplified  purpose,  the 
hot  and  cold  sources  are  assumed  to  be  infinite  constant  heat  res¬ 
ervoirs  in  the  present  model,  hence,  the  third  kind  boundary  con¬ 
dition  is  substituted  by  the  first  kind  one.  In  addition,  our 
preliminary  tests  showed  that  there  is  no  qualitative  difference  for 
dynamic  response  characteristics  of  the  TEG  although  quantitative 
difference  exists  when  this  simplification  is  adopted.  Adiabatic 
boundary  conditions  are  applied  to  the  outside  surfaces  of  the  TEG 
with  continuity  assumption  for  temperature  and  heat  flux  on  the 
interfaces  between  different  materials. 

For  the  electric  boundary  conditions  of  the  TEG,  assumptions 
of  constant  current  and  zero  electric  potential  are  adopted  for 
the  inlet  and  the  outlet,  respectively,  as  shown  in  Fig.  1.  There  is 
no  current  flowing  out  on  the  other  surfaces,  thus,  J  Ji  =  0. 
Besides,  initial  temperature  of  the  TEG  element  is  equal  to  the 
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Table  1 

Properties  of  the  materials  used  in  the  present  work  [28]. 


Materials 

k  (W  K1  rrr1) 

1/(7  (S  m  ^ 

a  (V  K-1) 

n-type  (Bi2Te3) 

1.6 

7.69  x  104 

-2.28  x  10"4 

p-type  (Sb2Te3) 

2.1 

9.62  x  104 

1.71  x  10"4 

Connector  (Cu) 

350 

5.90  x  108 

- 

Ceramic  plate  (Si) 

130 

- 

- 

environmental  temperature  Too  =  300  K,  and  initial  electric  po¬ 
tential  is  zero. 

3.  Simulation  cases 

The  geometry  of  the  TEG  adopted  for  simulations  has 
Ho  =  0.9  mm,  Hi  =  0.2  mm,  H2  =  1.9  mm,  L\  =  0.5  mm,  and 
L2  =  1.0  mm.  Sb2Te3  and  Bi2Te3  are  chosen  for  the  p-type  and  n-type 
semiconductors  due  to  their  excellent  thermoelectric  properties 
from  room  temperature  to  500  I<  [27].  The  connector  and  ceramic 
plate  are  made  of  copper  and  Si,  respectively.  All  of  the  material 
properties  are  listed  in  Table  1. 

The  transient  behaviors  of  the  TEG  are  mainly  trigged  by 
changes  of  the  working  environment  and  consumer’s  consumption 
demand  for  electric  power.  For  comparison  purpose,  a  base-line 
case  with  hot  end  temperature  Th  =  450  K,  cold  end  temperature 
Tc  =  300  K,  and  load  current  /  =  0.5  A  is  chosen  and  analyzed  firstly. 
The  step  increase  or  decrease  and  the  linear  increase  or  decrease  is 
designed  to  simulate  the  changes  of  Th,  Tc  and  /,  respectively.  Thus, 
twelve  cases  are  simulated  as  listed  in  Table  2,  in  which  only  one 
parameter  is  changed  for  each  case  with  the  other  two  parameters 
as  the  same  as  the  base-line  case. 

4.  Model  validation 

Due  to  the  lack  of  experimental  data  of  the  transient  response 
performance  of  the  TEG  in  the  open  literature,  the  model  is  vali¬ 
dated  only  for  the  steady-state  operation.  Qh,  l  P,  and  rj ,  were 
compared  with  one-dimensional  steady-state  analytical  solution  of 
Angrist  et  al.  [29]  and  numerical  prediction  of  Chen  et  al.  [8].  In  the 
analytical  solution,  the  thermal  conductivity,  electric  conductivity, 
and  Seebeck  coefficient  of  semiconductors  were  all  assumed  to  be 
temperature-independent,  and  adiabatic  boundary  conditions 
were  assumed  on  the  side  surfaces  of  TEG  element.  Our  simulation 
is  performed  for  a  TEG  device  with  127  elements,  each  element 
has  the  geometry  of  Hi  =  0.2  mm,  H2  =  1.6  mm,  Li  =  0.5  mm, 
and  L2  =  1.4  mm,  as  well  as  constant  material  properties  of 
«p  =  -an  =  226.8  x  10”6  V  I<-1,  kp  =  kn  =  1.52  W  m”1  K_1  and 


Table  2 

Cases  simulated  in  the  present  work. 


Case 

Parameters  Value  before 

changing 

Variational 

form 

Value  after  Unit 

changing 

1 

Th 

450 

500 

K 

2 

Tc 

300 

_r 

320 

K 

3 

I 

0.5 

0.7 

A 

4 

Th 

450 

500 

K 

5 

Tc 

300 

— / 

320 

K 

6 

I 

0.5 

0.7 

A 

7 

Th 

450 

400 

K 

8 

Tc 

300 

“L 

280 

K 

9 

I 

0.5 

0.3 

A 

10 

Th 

450 

400 

K 

11 

Tc 

300 

280 

K 

12 

I 

0.5 

0.3 

A 

(7p  =  {7n  =  1.447  x  1CT5  Q  m.  In  addition,  the  heat  loss  to  the  ambient 
is  ignored  for  comparison  with  the  theoretical  solution.  The  same 
geometry  and  material  properties  were  also  used  in  Chen  et  al.’s 
simulation  [8]. 

The  simulation  was  conducted  for  three  sets  of  grid  systems 
(coarse  grid,  medium  grid,  and  refined  grid).  The  present  and  Chen 
et  al.’s  [8]  numerical  predictions  and  the  theoretical  evaluation  [29] 
for  Qh,  /,  P,  and  rj  are  shown  in  Table  3.  Present  prediction  has  better 
agreement  with  the  theoretical  evaluation  than  Chen  et  al.’s.  The 
maximum  deviation  between  the  present  prediction  with  the 
coarse  grid  and  the  theoretical  evaluation  is  3.5%,  2.5%  for  the 
medium  grid,  and  2.0%  for  the  refined  grid.  Restated  that  the 
theoretical  solution  [29]  was  driven  based  on  the  assumptions  of 
steady-state  operation,  one-dimensional  temperature  distribution, 
uniform  current  density,  constant  material  properties,  and  no  heat 
loss  to  the  ambient,  therefore,  its  applicability  will  be  limited  when 
the  operating  condition  of  the  TEG  deviates  from  these  assump¬ 
tions.  However,  the  present  model  is  more  general  with  three- 
dimensional  distributions  of  temperature  and  current  density, 
which  can  be  used  to  study  the  transient  or  steady  performance  of 
the  TEG  at  more  complicated  operating  conditions  such  as  high 
temperature  difference  and/or  large  heat  loss  to  the  ambient. 

The  grid  independence  for  the  TEG  geometry  used  in  the  case 
studies  was  examined  in  preliminary  test  runs.  The  independence 
of  the  time  step  size,  At,  was  also  examined  as  shown  in  Fig.  2  with 
At  then  taken  as  0.001  s. 

5.  Results  and  discussion 

5.1.  The  effect  of  hot  end  temperature 

5.1.1.  Step  and  linear  decreases 

Fig.  3  shows  variations  of  P  with  time  for  the  step  and  linear 
decreases  in  Th.  For  the  step  variation,  Th  decreases  sharply  from 
450  K  to  400 1<  at  0.01  s,  however,  P  does  not  change  synchronously 
which  reaches  the  new  steady-state  value  of  0.0083  W  at 
t  =  0.052  s.  For  the  linear  variation,  Th  decreases  from  450  K  at 
t  =  0.01  s  to  400  K  at  t  =  0.06  s  with  reduction  rate  of  1  K  m  s_1,  and 
the  new  steady  state  is  observed  at  t  =  0.086  s  with  the  same  P 
value  of  0.0083  W.  The  results  indicate  that  there  exists  a  temporal 
hysteresis  in  the  response  of  P  to  Th.  The  hysteresis  is  attributed  to 
the  delay  of  thermal  diffusion  from  the  ceramic  plate  to  the  hot 
junction.  The  response  time  is  defined  as  the  time  required  to  reach 
a  new  steady  state  after  the  outside  environment  (Th,  Tc  or  /)  stops 
changing.  The  response  time  of  P  is  0.042  s  (=0.052  s  -  0.01  s)  for 
the  step  decrease  and  0.026  s  (=0.086  s  -  0.06  s)  for  the  linear 
decrease,  which  indicates  that  the  response  time  shortens  once  the 
reduction  rate  of  Th  is  reduced. 

Fig.  4  shows  the  dynamic  variations  of  Qh  for  the  step  and  linear 
decreases  in  Th.  For  the  step  variation,  Th  decreases  from  450  K  to 
400  K  at  t  =  0.01  s,  which  causes  a  sharp  reduction  in  Qh  from 
0.3721  W  to  -68.4293  W  at  t  =  0.011  s,  then  Qh  gradually  increases 
and  again  reaches  zero  at  t  =  0.051  s,  finally,  the  new  steady  state 

Table  3 

Comparisons  of  numerical  predictions  and  analytical  solutions  (Th  =  423  I<  and 
Tc  =  303  K). 

Quantity  Analytical  Chen  et  al.’s  Grid  systeml  Grid  system2  Grid  system3 
results  [29]  results  [8]  (75081  grids)  (241101  (557289 

grids)  grids) 


Qh  (W) 

81.3 

81.8 

81.2 

81.2 

81.2 

P(W) 

3.98 

3.62 

3.84 

3.88 

3.90 

1(A) 

1.08 

1.032 

1.08 

1.08 

1.08 

7] 

4.89 

4.43 

4.73 

4.77 

4.80 
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Fig.  2.  Dynamic  variations  of  Qh  and  Qc  with  various  time  step  lengths  for  increased 
load  current. 


corresponding  to  Th  =  400  K  is  reconstructed  at  t  =  0.095  s.  For  the 
linear  variation,  Qh  firstly  decreases  from  the  initial  value  of 
0.3721  W  to  the  minimum  value  of  -5.8711  W  (Qh.min)  at  t  =  0.06  s, 
then  increases  to  zero  at  t  =  0.085  s  and  finally  reaches  the  steady 
state  at  t  =  0.127  s.  It  is  noted  that  a  negative  value  of  Qh  occurs 
during  0.01-0.052  s  for  the  step  variation  and  0.01-0.085  s  for  the 
linear  variation,  which  means  that  the  hot  end  of  the  TEG  does  not 
absorb  but  liberate  the  heat  to  the  ambient,  however,  the  TEG  still 
generates  electricity.  This  abnormal  phenomenon  can  be  explained 
by  the  internal  temperature  distribution  of  the  TEG  as  shown  in 
Fig.  5. 

Fig.  5  shows  the  temperature  distributions  along  the  TEG  height 
direction  at  x  =  0.75  mm  and  at  different  times  for  the  step  and 
linear  decreases  in  Th.  For  the  step  variation,  Th  is  450  I<  and  the 
temperature  at  the  top  of  p-type  semiconductor,  TPit0p,  is  448.0  K  at 
t  =  0.01  s,  with  only  2  I<  temperature  difference  due  to  high  thermal 
conductivities  of  the  ceramic  plate  and  metal  connector.  At 
t  =  0.011  s,  Th  decreases  sharply  to  400  K,  however,  TPit0p  is  only 
decreases  to  443.6  I<  due  to  the  delay  of  the  thermal  diffusion  from 
the  ceramic  plate  to  the  hot  junction.  TP)t0p  gradually  decreases  at 
0.011  s  <  t  <  0.052  s,  however,  it  always  higher  than  Th,  thus,  a 
positive  temperature  gradient  from  the  hot  junction  to  the  ceramic 
plate  is  formed,  which  causes  a  heat  release  of  the  TEG  to  the 
ambient.  The  maximum  temperature  gradient  occurs  at  t  =  0.011  s 


Fig.  3.  Dynamic  output  power  variations  for  the  linear  and  step  decreases  in  Th. 


for  the  step  variation,  so  that  the  maximum  heat  release  is  observed 
at  the  same  moment  (Fig.  4).  At  t  =  0.052  s,  TP(t0p  decreases  to 
399.2  K,  which  is  lower  than  Th,  and  hence  the  TEG  starts  to  absorb 
the  heat  from  the  ambient  again.  P  is  mainly  determined  by 


Fig.  5.  The  temperature  distributions  along  the  TEG  height  direction  at  different  times: 
(a)  step  decrease:  (b)  linear  decrease. 
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temperature  distribution  of  p-type  and  n-type  semiconductors. 
Fig.  5a  shows  that  the  temperature  at  the  bottom  of  the  p-type 
semiconductor,  TPib0ttom,  is  almost  kept  at  300  K,  while  TP)t0p 
continuously  decreases  from  448.0  I<  (t  =  0.01  s)  to  398.6  K 
(t  =  0.095  s)  throughout  the  whole  transient  process,  which  causes 
a  continuous  reduction  in  the  temperature  difference  between 
Tp, bottom  and  TP)top  and  hence  a  continuous  reduction  in  P  of  the  TEG 
(Fig.  3).  For  the  linear  variation,  Th  decreases  from  450  K  to  400  K 
with  a  finite  reduction  rate  of  1  I(  m  s-1,  the  response  of  TP)top  is 
slower  than  the  step  variation.  For  example,  TP)t0p  decreases  to 
409.3  I<  at  t  =  0.025  s  for  the  step  variation  (Fig.  5a),  however,  it  only 
decreases  to  432.4  I<  for  the  linear  variation  (Fig.  5b).  The  slower 
response  of  TP)t0p  is  the  reason  why  the  times,  at  which  Qh.min, 
Qh  =  0  and  the  new  steady  state  are  reached,  are  delayed  for  the 
linear  variation  as  compared  to  the  step  variation. 

The  dynamic  variations  of  77  for  the  step  and  linear  decreases  in 
Th  are  shown  in  Fig.  6.  According  to  Eq.  (7),  77  is  defined  as  the  ratio 
of  P  to  Qh-  At  t  <  0.01  s,  the  TEG  operates  at  steady  state  with 
Th  =  450  I<  and  Tc  =  300  K,  and  77  is  4.85%.  At  0.01  s  <  t  <  0.051  s  for 
the  step  variation  and  at  0.01  s  <  t  <  0.085  s  for  the  linear  variation, 
the  TEG  liberates  the  heat  to  the  ambient,  which  implies  that  the 
output  power  is  generated  without  the  outside  heat  input,  so 
that  an  infinite  77  occurs.  At  t  =  0.052  s  for  the  step  variation  and 
at  t  =  0.086  s  for  the  linear  variation,  the  TEG  starts  to  absorb 
the  heat  from  the  ambient  again,  thus,  a  finite  77  =  31.69%  for 
the  step  variation  and  77  =  39.72%  for  the  linear  variation  occurs, 
then  77  decreases  due  to  the  reduced  temperature  difference 
(=7p,top  -  Tp, bottom)  and  P  until  the  TEG  reaches  the  new  steady  state 
with  77  =  3.12%. 


Fig.  7.  Dynamic  variations  of  the  power  and  conversion  efficiency  for  the  step  and 
linear  increases  in  Th. 


Fig.  7  also  shows  that  there  exists  an  undershoot  of  77  for  both 
step  and  linear  increases  despite  the  continuously  increased  P.  The 
undershoot  occurs  at  0.01  s  <t<  0.011  s  for  the  step  increase  and  at 
0.01  s  <  t  <  0.06  s  for  the  linear  increase.  The  undershoot  is  more 
intense  for  the  step  increase,  77  almost  reduces  to  zero  at  t  =  0.011  s, 
which  is  much  lower  than  its  initial  value  of  0.0180  W.  This  phe¬ 
nomenon  can  be  attributed  to  the  lag  of  TP)t0p  ( ~  Tn  top)  variation 
behind  Th  due  to  finite  thermal  diffusion  rate  from  the  ceramic  plate 
to  the  hot  junction  of  the  TEG. 


5.1.2.  Step  and  linear  increases 

Dynamic  variations  of  the  P  and  77  for  the  step  and  linear  in¬ 
creases  in  Th  are  shown  in  Fig.  7.  For  the  step  variation,  Th  increases 
sharply  from  450  K  to  500  I<  at  t  >  0.01  s,  however,  for  the  linear 
variation,  Th  increases  from  450  I<  to  500  I<  at  0.01  s  <  t  <  0.06  s 
with  an  increase  rate  of  1  I<  m  s-1.  Because  the  increase  in  Th  ele¬ 
vates  the  temperature  difference  between  the  hot  and  cold  junc¬ 
tions,  Tp  top  -  Tpbottom  ( ~  Tf^top  -  Tnbottom)>  ^nd  the  temperature 
difference  becomes  larger  and  larger  with  the  time  increased,  P 
gradually  rises  until  the  TEG  reaches  the  new  steady  state.  The 
temporal  hysteresis  in  the  response  of  P  to  Th  is  also  observed  for 
the  present  two  cases  due  to  the  delay  of  the  thermal  diffusion  from 
the  ceramic  plate  to  the  hot  junction. 


in  rh- 


5.2.  The  effect  of  cold  end  temperature 

Step  and  linear  variations  in  Tc  are  designed,  respectively,  to 
analyze  the  effect  of  Tc  on  the  dynamic  characteristics  of  the  TEG. 
The  TEG  operates  at  steady  state  with  Th  =  450  K  and  Tc  =  300  K  at 
t  <  0.01  s,  then  Tc  starts  to  increase  or  decrease  from  300  K  with 
AT  =  20  K  at  t  >  0.01  s.  For  the  linear  variation,  the  increase  or 
decrease  rate  of  Tc  is  taken  as  0.4  K  m  s-1.  Fig.  8  shows  that  the 
power  curves  for  the  increase  and  decrease  cases  have  symmetry. 
For  example,  P  is  increased  from  its  initial  value  of  0.0180  W— 
0.0219  W  with  AP  =  0.039  W  at  0.038  s  for  the  step  decrease,  while 
is  decreased  to  0.0141  W  with  the  same  AP  =  0.039  W  at  the  same 
time  for  the  step  increase.  Reduction  of  P  due  to  increasing  Tc  can  be 
attributed  to  reduced  temperature  difference  between  the  hot  and 
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cold  ends  of  the  TEG,  and  vice  versa.  The  hysteresis  of  response  of  P 
to  Tc  is  once  again  observed  for  all  the  four  cold  temperature  vari¬ 
ations.  The  TEG  reaches  the  new  steady  state  at  0.038  s  with  the 
response  time  of  0.028  s  for  the  two  step  variations,  while  the 
response  time  shortens  to  0.013  s  (=0.073  -  0.06)  for  the  two  linear 
variations. 

Dynamic  variations  of  tj  for  the  step  and  linear  variations  in  Tc  is 
shown  in  Fig.  9.  Note  that  overshoot  of  rj  is  observed  for  the  step  and 
linear  decreases  in  Tc,  while  undershoot  for  the  step  and  linear 
increases.  Due  to  the  delay  of  thermal  diffusion,  Th  has  the  slowest 
response  to  Tc,  and  then  followed  by  Tpit0p  (~Tn,t0p)  and  TPibottom 
( ~  Tn> bottom)-  Thus,  for  the  decrease  cases  in  Tc,  the  increase  rate  of 
Qh  is  lower  than  P,  which  causes  the  overshoot  of  77.  The  response  of 
TPi bottom  (~Tn, bottom)  for  the  step  decrease  is  faster  and  stronger 
than  that  for  the  linear  decrease,  and  hence  the  overshoot  of  77 
occurs  earlier  with  the  higher  peak  value.  The  similar  reason  can  be 
used  to  explain  the  undershoot  of  77  for  the  increase  cases  in  Tc. 

5.3.  The  effect  of  load  current 

The  simulations  are  performed  with  fixed  Th  =  450  K  and 
Tc  =  300  K.  The  initial  load  current,  /,  is  0.5  A  at  t  <  0.01  s,  it  varies  to 
0.7  A  at  t  >  0.01  s  for  the  step  and  linear  increases,  and  to  0.3  A  for 


the  step  and  linear  decreases.  The  increase  or  decrease  rate  of  I  is 
4  A  s-1  for  the  linear  cases.  The  dynamic  responses  of  Qh  for  the  step 
and  linear  variations  in  I  are  shown  in  Fig.  10.  The  steady  state  is 
reached  at  0.039  s  for  the  step  increase,  at  0.043  s  for  the  step 
decrease,  at  0.075  s  for  the  linear  increase,  and  at  0.080  s  for  the 
linear  decrease.  These  results  indicate  that  Qh  lags  behind  I.  The  lag 
behavior  can  be  explained  as  follow.  Because  Joule  heat  and 
Thomson  heat  are  volumetric  heat  sources  and  are  proportional  to 
I2  and  /,  although  the  internal  heat  sources  response  immediately 
when  I  is  varied,  the  temperature  distribution  can  not  response 
synchronously  due  to  the  finite  thermal  diffusion  rate. 

Dynamic  variations  of  P  for  the  step  and  linear  cases  are  shown 
in  Fig.  11.  The  response  of  P  is  almost  synchronous  with  /,  which 
indicates  that  the  electric  response  of  the  TEG  is  much  faster  than 
the  thermal  response.  The  steady-state  I—P  curve  with  fixed 
Th  =  450  K  and  Tc  =  300  K  is  shown  in  Fig.  12.  If  the  P—t  curve 
(0.01  s  <t<  0.06  s)  in  Fig.  11  is  converted  into  the  I—P  curve,  the  I—P 
curve  almost  coincide  to  the  steady-state  one  in  the  same  range  of 
load  current  (0.5  A  <  I  <  0.7  A)  as  shown  in  Fig.  12,  which  further 
confirms  that  the  electric  response  time  of  the  TEG  can  be  almost 
neglected.  Dynamic  variations  of  conversion  efficiency  for  the  step 
and  linear  variations  in  I  are  shown  in  Fig.  13.  Temporal  hysteresis  in 
the  response  of  77  to  /  for  the  four  cases,  overshoot  for  the  two 


1(A) 


Fig.  12.  Comparison  between  the  steady-state  and  dynamic  I-P  curves  with  fixed 
Th  -  450  K  and  Tc  =  300  K. 
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in  /. 

increase  cases,  and  undershoot  for  the  two  decrease  cases  are 
observed,  respectively.  These  phenomena  can  be  attributed  to  the 
synchronous  variation  of  P  and  the  lagging  of  Qh  when  I  is  varied. 

6.  Conclusions 

This  paper  developed  a  complete,  three-dimensional  and  tran¬ 
sient  TEG  model,  the  model  solved  the  temperature  and  electric 
potential  equations  and  accounted  for  all  thermoelectric  effects. 
The  model  was  used  to  investigate  the  dynamic  characteristics  of 
the  TEG  when  the  application  environments  are  changed.  Ther¬ 
moelectric  materials  Bi2Te3  and  Sb2Te3  were  selected  as  n-type  and 
p-type  semiconductors.  Four  kinds  of  typical  cases,  the  step  in¬ 
crease  and  decrease  and  the  linear  increase  and  decrease  in  the  hot 
end  temperature,  the  cold  end  temperature,  and  the  load  current, 
were  designed  to  investigate  their  effects  on  the  dynamic  charac¬ 
teristics  of  the  TEG.  The  main  conclusions  are  as  follows: 

1 )  For  the  variation  of  the  hot  end  temperature,  increasing  the  hot 
end  temperature  elevates  the  output  power  and  the  heat 
amount  adsorbed  from  the  ambient,  however,  the  growth  rate 
of  the  heat  amount  is  higher  than  that  of  the  output  power, 
which  causes  an  undershoot  of  the  conversion  efficiency. 
Decreasing  the  hot  end  temperature  reduces  the  output  power 
and  produces  a  positive  temperature  gradient  from  the  hot 
junction  to  the  ceramic  plate,  which  causes  a  heat  release  of  the 
TEG  to  the  ambient  and  infinite  conversion  efficiency.  The 
temporal  hysteresis  in  the  response  of  the  output  power  to  the 
hot  end  temperature  is  observed  for  all  the  four  cases  due  to  the 
delay  of  the  thermal  diffusion  from  the  ceramic  plate  to  the  hot 
junction. 

2)  When  the  cold  end  temperature  is  varied,  increasing  the  cold 
end  temperature  reduces  the  output  power,  and  vice  versa.  The 
response  hysteresis  of  the  output  power  to  the  cold  end 


temperature  is  observed  and  the  response  time  for  the  linear 
variation  is  smaller  than  the  step  variation.  The  conversion 
efficiency  displays  overshoot  for  the  decrease  cases  and  un¬ 
dershoot  for  the  increase  cases. 

3)  When  load  current  is  varied,  the  response  of  the  output  power 
is  almost  synchronous  with  the  load  current,  however,  the 
temperature  distribution  within  the  TEG  and  the  correspond¬ 
ing  heat  adsorbed  from  the  ambient  can  not  response  syn¬ 
chronously  due  to  finite  thermal  diffusion  rate,  which  causes 
not  only  the  temporal  hysteresis  in  the  response  of  conversion 
efficiency  to  the  load  current,  but  also  overshoot  or  undershoot 
of  the  conversion  efficiency. 
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